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' We check the impact of the factorization scheme employed in the calculation of the 

, heavy-quark deep-inelastic scattering (DIS) electro-production on the PDFs determined 

1 I . in the NNLO QCD analysis of the world inclusive neutral-current DIS data combined 
i-^ ' with the ones on the neutrino-nucleon DIS di-muon production and the fixed-target 

Drell-Yan process. The charm-quark DIS contribution is calculated in the general- 

^ • mass variable-flavor-number (GMVFN) scheme: At asymptotically large values of the 

' momentum transfer Q it is given by the zero-mass 4-flavor scheme and at the value 

. of Q equal to the charm-quark mass it is smoothly matched with the 3-flavor scheme 

' using the Buza-Matiounine-Smith-van Neerven prescription. The PDFs obtained in 

I I this variant of the fit are very similar to the ones obtained in the fit with a 3-fiavor 

^ ■ scheme employed. Our 5-flavor PDFs derived from the 3-flavor ones using the NNLO 

QQ ' matching conditions are used to calculate the rates of W^/Z and production at the 

I Tevatron coUider and the LHC at NNLO. 

^■f^ . 1. A clear theoretical understanding of the heavy-quark electro-production is necessary 

I for the phenomenology of lepton-nucleon DIS at HERA since it contributes up to 25% to the 

«) i DIS inclusive cross section at small a; [2l [3] . This is particularly important for the validation 

' of the parton distribution functions (PDFs) in the kinematic region relevant for the foreseen 

. experiments at the LHC since, despite the impact of collider data on the PDFs, DIS remains 

^ ' an unique source of information about PDFs at small x. Two complementary approaches 

were suggested to calculate the DIS heavy-quark contribution. One of them is based on 
the factorization scheme with three light quarks in the initial state. In the 3-flavor scheme 
the heavy-quark production is calculated in fixed order of perturbative QCD with complete 
corrections up to the 0{al) and partial corrections up to 0{al) available [H [HI [H [7] . The 
logarithmic terms in the QCD corrections rise in the asymptotic region of the momentum 
transferred Q 3> rrih-, where uih is the heavy quark mass. Therefore they must be re-summed 
for such kinematics [8] . This leads to the concept of the variable-flavor scheme ( VFN) with 
the heavy quarks considered in the same way as the light-flavor ones. In the VFN scheme 
large-log re-summation suggested in Ref. [5] is performed through the QCD evolution of the 
heavy-quark PDFs [5] . The VFN scheme is valid only at asymptotically large values of Q 
and cannot be routinely extrapolated to the low-Q region. For realistic kinematics it has to 
be extended to the case of a general-mass VFN (GMVFN) scheme, which, at small Q, must 
give a smooth matching with the 3-flavor scheme. Different variants of the GMVFN scheme 
are used in the context of the global flt of PDFs based on small- a; DIS data [TUj. At the 
matching point Q = rrih the structure function F2 calculated within the GMVFN schemes 
of Ref. |10| is continuous. However these schemes do not guarantee a smooth transition of 
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Figure 1: The values of F2{x,Q^) calculated in the different schemes compared to the data by 
HI [2] (squares) and ZEUS (circles) collaborations (dotted-dashes: the 3-flavor scheme; solid 
lines: the BMSN prescription of the GMVFN scheme, dashes: the massless 4- flavor scheme). 

F2 between the 3-flavor and the VFN schemes. In our study we consider a variant of the 
GMVFN scheme based on the Buza-Matiounine-Smith-van Ncerven (BMSN) prescription 
of Ref. [TT], which indeed provides a smooth behavior of Fj* at Q = rrih, and perform a 
phenomenological comparison of this prescription with the 3-flavor scheme in the PDFs fit. 

2. The conclusion about the impact of the scheme evidently depends on the kinematic 
coverage and precision of the data used in the fit. Figure [1] shows a comparison of the 
charm structure function F2 calculated in the 3-fiavor scheme and the BMSN prescription 
of the GMVFN scheme with the use of the MRST2001 PDFs [Hj to the representative set 
of the recent HERA collider data by HI and ZEUS [HIS]. The available data on F| are 
not sensitive to the difference between the 3-flavor and the GMVFN scheme in the BMSN 
prescription. The best sensitivity of the data to the scheme choice appears at x ~ 0.003, 
however even in this case the errors in the data are quite big as compared to the scheme 
differences. This docs not look like a specific feature of the BMSN prescription since even at 
maximal values of Q reached by the existing experiments the difference between the massless 
4-fiavor scheme and the 3-fiavor scheme is smaller than the data uncertainties. In particular 
this happens because the NLO corrections of Ref. [5| to the heavy-quark electro-production 
greatly reduce the need of large- log re-summation [T^. With the NNLO corrections of 
Ref. [7] this need should be further reduced. At small Q the difference between the 3-flavor 
and BMSN scheme is suppressed, similarly to any other GMVFN scheme, which provides 
smooth matching with the former. The discrepancy between theory and data at small Q and 
X observed in Figure [1] cannot be resolved by a scheme change. It rather can be cured by the 
partial NNLO corrections of Ref. [M] due to the soft gluon re-summation near threshold. 

3. The inclusive DIS data are more sensitive to the choice between the VFN and 3-flavor 
scheme because of higher accuracy. However, the PDFs obtained from the fits based on 
these two variants are still similar. We check this point for the example of a combined fit 
of PDFs to the inclusive DIS data supplemented by the fixed-target Drell-Yan data and 
the DIS di-muon production data, which improves the fiavor separation of the quark PDFs. 
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Figure 2: The gluon (left panel) and sea (right panel) distributions at the factorization scale of 
fi = 3 GeV obtained in two variants of the fit. Solid lines; 3-flavor scheme, dashed lines ; GMVFN 
scheme in the BMSN prescription. 

We consider two variants of the fit, with F| calculated in the 3-flavor scheme or in the 
BMSN prescription. The NNLO corrections to the PDF evolution [TS], the light-parton 
coefficient functions [Hj, and the NLO corrections of Ref. [S] to the heavy-quark electro- 
production coefficient functions are taken into account. The non-QCD corrections applied 
and other details of the fit can be found elsewhere [T7|. The gluon and sea distributions 
obtained in the two variants of the fit are compared in Figure [H The difference between 
the two variants of the fit is marginal and is in agreement with the comparison given in 
Figure [1] In the case of Thorne-Roberts (TR) prescription of Ref. [18] employed for the 
GMVFN scheme the impact of the scheme choice is somewhat bigger [19]. That might be 
related to a kink in at Q = rric for small x. For the ACOT prescription jp] and the 
ACOT(x) modification of this prescription [20] employed in the NNLO analysis of Ref. [21] 
this kink is more pronounced than for the case of the TR prescription. As a result, the heavy- 
quark contribution to the inclusive structure function is overestimated and the light-parton 
distributions arc underestimated, correspondingly. 

4. Contrary to the case of DIS for many processes of interest at the energies of the 
Tevatron collider and the LHC the use of a VFN scheme is justified since the typical fac- 
torization scale is much bigger than the c— and &— quark masses. Moreover, very often such 
an approach is very efficient since it allows to simplify matrix element calculations. For 
the purposes of collider phenomenology applications we generate the NNLO 4-flavor PDFs 
using the matching conditions of Ref. [11] with the 3-flavor PDFs obtained in our fit as 
input. Similarly, we use these 4-flavor PDFs to produce the 5-flavor ones. Figure [3] shows 
a comparison of our 4-flavor PDFs at the scale of ^ = 2 GeV with the NNLO MSTW08 
PDFs of Ref. [H]. At this scale the sea and gluon distributions by MSTW are smaller 



"At the factorization scales close to rric the 4-fiavor PDFs arc inapplicable. However, we use this set for 
the purpose of comparison since the 3-flavor PDFs are unavailable for the MSTW08 case. 
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Figure 3: The la band for the representative set of our PDFs (shaded area) and the central values 
of other NNLO PDFs sets (dashes). Left panel: the 3-flavor u-, d-quarks and giuons at the scale of 
fi = 2 GeV compared to ones of the JR set; central panel: the 4-flavor u-, d-quarks and giuons at 
the scale oi n = 2 GeV compared to ones of the MSTW08 set; right panel: the 5-flavor s-, c-, and 
b-quarks at the scale of /i = 100 GeV compared to ones of the MSTW08 set. 



than ours, moreover the MSTW gluon distributions are negative at a; < 0.0005. The dif- 
ference in the gluon distribution of these two sets is sensitive to the recent measurement 
of the structure function Fl by the HI collaboration [23]. The HI data are in good agree- 
ment to the predictions of based on our PDFs, however, then go slightly above the 
MSTW08 ones and in that way prefer bigger giuons at small x. The agreement between our 
PDFs and the JR PDFs of Ref. [22] obtained in the fit performed in the 3-flavor scheme 

is much better (see Figure [3]). In this 
case the gluon distributions differ only 
at a; 10"^, in a region not being 
constrained by data. Our 5-flavor c— 
and 6— quark distributions arc larger 
than the MSTW ones. This is related 
to the difference in the gluon distribu- 
Table 1: The integral VV'^/Z and ti production cross tions. The uncertainties in our PDFs 
sections at the energies of Tevatron and LHC with the given in Figure [3] are propagated from 
la uncertainties estimated from the fit resuhs in the ^^^g uncertainties in experimental data 
present analysis. ^-^^ account of their correlations. We 

also include the effect of the c— and b— 
quark mass variation by ±0.1 GeV and 
±0.5 GeV, respectively, into the PDF uncertainties. For the case of c— and 6— quark PDFs 
the uncertainties in the masses are dominant sources. The value of the strong coupling 
as{Mz) = 0.1135 ± 0.0014 is obtained in our fit, in good agreement with the results of 
Ref. [12 and earher results of Refs. fM[\Tl\. 

5. The rates of the /Z and tt production at the Tevatron collider and the LHG 
obtained with our 5-flavor PDFs are given in Table 1. The W^/Z rates are calculated 
with the QGD corrections up to NNLO [25]. In the case of the tt production only partial 



V^s (TeV) 


a{W^) [nb] 


a{Z) [nb] 


a{tt) [pb] 


1.96 (pp) 


26.2 ±0.3 


7.73 ± 0.08 


6.91 ±0.17 


7 (pp) 


98.8 ± 1.5 


28.6 ±0.5 


131.3 ± 7.5 


10 (pp) 


145.6 ±2.4 


42.7 ±0.7 


342 ± 15 


14 (pp) 


207.4 ±3.7 


61.4 ± 1.1 


780 ± 28 
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NNLO corrections of Ref. [55] stemming from the threshold ghion re-summation are taken 
into account. The uncertainties in the /Z rates due to the PDFs vary from 1% at the 
Tevatron energy to 2% at the maximal LHC energy. Such a level of accuracy allows to 
use these processes as a luminosity monitor. Our estimate of the tt production rate for 
Tevatron is in agreement to the estimate obtained with the MSTW08 PDFs, while at the 
LHC energies the latter are bigger than ours. As in the case of Fl this difference is related 
to the difference in the gluon distributions. Therefore tt production at the LHC can be used 
to discriminate between these two PDFs sets. 
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